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Multifragmentation reactions without large compression in the ini- 
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Q . tial state (proton-induced reactions, reverse-kinematics, projectile frag- 

mentation) are examined, and it is shown that the high temperatures 
obtained from fragment kinetic energy spectra and lower temperatures 

Q\ . obtained from observables such as level population or isotope ratios 

can be understood in a common framework. 



The phenomenon of nuclear niultifragmentation, the decay of highly ex- 
cited nuclear matter into (nucleons and) several fragments of size A^ > 5, 
has been observed in proton-induced and heavy ion reactions |I|, |^. Very 
recently, similar multifragmentation events have also been observed in the 
disintegration of Cgo-fullerenes ('buckyballs') after bombardment with high- 
energy heavy ions . It is clear by now that there is probably no one single 
mechanism responsible for the multitude of different fragmentation scenar- 
ios which range from sequential binary emission processes at low excitation 
energies [Q over statistical multifragmentation f^ to explosive vaporization. 
The probably most exciting possibility is the occurrence of a 2'^'^-order phase 
transition at the critical point of nuclear matter @, |^, §]. 

The majority of theoretical investigations has been aimed at understand- 
ing the fragment mass distributions. However, there is considerable interest 
in fragment kinetic energy spectra, and the present letter is also concerned 
with this topic. In high-energy proton-induced ^ and in relativistic projec- 
tile fragmentation ^ reactions one observes Boltzmann-like kinetic energy 
spectra of intermediate mass fragments with temperatures, for which the frag- 
ment mass dependence has been parameterized [§] as Tf = To {A,- — Af)/^^ 
where A,- is the mass of the recoil residue. The temperature To typically 
has a value of ~ 15 MeV, much larger that nuclear temperatures extracted 
from isotope ratios or level population ratios |]10|], and also much larger than 



typical temperatures (~ 5 — 8 MeV) used in statistical models used to 
reproduced the experimental mass yield curves. 

Here it is argued that the comparatively high value of Tf is a conse- 
quence of the addition of the Fermi momenta of the individual nucleons in 
the fragment and thus a consequence of the Fermi-Dirac nature of nucle- 
ons. The single-particle model is employed. It previously was successfully 



used ||TT], |T2|, |I^ to explain the observed dispersion in fragment transverse 
momentum spectra generated in projectile fragmentation. 

One starts with an ensemble of A nucleons at (internal) temperature 
Tin. In the single-particle picture their momentum distribution (in the non- 
relativistic limit) is given by 

p(p) = (1 + exp[(pV2m - ^)/TiJ)-i , (1) 

where /i is the chemical potential. Adding up the individual momentum 
vectors of all A^ nucleons in a fragment one then arrives at the momentum 



distribution of the fragment 

p(Pf) = / Eid'p^ p(p.)}5'(Pf - E P.) • (2) 

This addition procedure is of course nothing else but a random walk in 
momentum space and thus in the class of problems first posed by Pearson 
[Q. Since the single-particle probability distribution has average and a 
fixed variance a^, the central limit theorem of Gauss applies in the present 
case, and - for sufficiently large number of steps, A^, in the random walk - 
the probability distribution can be written as 

Using the non-relativistic approximation E^ = P^/{2mi\fA{), we therefore 
find an exponentially falling kinetic energy spectrum for the fragments 

p(Ef) = ^=^^exp (-^) , (4) 



where in this picture the apparent temperature, Tf, is given by 

Tf = a'^/mN ■ (5) 

What exactly constitutes a sufficiently large number of steps is not a 
priory clear. However, in Fig. 1 I answer this question numerically by solving 
Eq. ^ with a dispersion of cr^ = {2/b)EpmM = 0.014 GeV^. (This choice 
will become more apparent below). For each fragment mass 5 x 10^ events 
were generated. Displayed are the resulting distributions (histograms) for 
Af = 2, 3, 5, 12, 100. One can clearly see that even a very small number of 
steps (between 5 and 12) are sufficient to approach the limit given by Gauss, 
Eq. I (sohd line). 

In the simple Fermi-gas model at temperature, the variance in the 
momentum distribution is given by o"^ = p'p/5, where pp is the Fermi mo- 
mentum. However, as pointed out by Goldhaber fl^ , there is the additional 
constraint that the momenta of all A nucleons combined add up to 0. This 
leads to 
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2 ^ ~ ^f / 2 
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(Pt) ■ (6) 



The term {A — Af)/{A — 1) is the recoil correction due to total momentum 
conservation. It explains the experimentally observed dependence of Tf on 
Af §. 

Inserting Eq. ^ into Eq. ^ leads to an apparent fragment kinetic temper- 
ature 

1 A-Af ^ A-Af2 

where (Ek) is the average kinetic energy per nucleon as obtained from the 
single particle model. In the finite temperature case considered here this 
average is 

^^'^ ^ Jo ^'l + exp[(e-^)/ri„] / Jo '^'l + exp[(6-^)/ri„] ' ^^^ 

In the limit Tin —>■ 0, we obtain /i —^ Ep and (1 + exp[(e — fi)/Tir])~'^ — *■ 
9{Ep — e), and the average of Eq. |^ becomes simply {Ek{Ti^ = 0)) = ^Ep, 
and consequently, the apparent fragment kinetic temperature at T[^ = will 
not be 0, but rather 

T,iT,^ = 0) = ^^lEp. (9) 

In the case Tin > 0, the single particle model is still applicable. Its 
assumptions are not dependent on any sudden approximation (as frequently 
applied in the study of projectile fragmentation), but also hold for a system 
in equilibrium at some finite temperature. This was already pointed out 
by Goldhaber |]T2[. In this case, the numerical solution of Eq. |^ inserted 
into Eq. ^ will yield the desired answer, the dependence of the apparent 
temperature as extracted from fragment kinetic energy spectra on the 'real' 
temperature of the system. For small temperatures. Tin, on can perform a 
Sommerfeld expansion of these equations around Tin = 0. This yields the 
approximation 

for the chemical potential and for the apparent fragment temperature the 
useful expression 

T,^T,^^^^^%(l^'4(M^Om'] , (11) 



which works to better than five percent up to Ti^/Ep = 0.3. Fig. 2 shows 
a comparison of the numerical solution (ignoring the recoil correction term 
{A — Af)/{A — 1)) (solid line) as well as the approximation of Eq. |TT[ 

Bertsch [jl3[ has pointed out that there is (for Tjn = 0) a suppression of 



the dispersion due to Pauli correlations, which can amount to about 30% in 
the projectile fragmentation of ^°Ca. This effect is not included in the simple 
estimate provided here. 

The effects of final state interaction, the emission of a nucleon from the 
excited pre-fragment, on the fragment kinetic energy distribution is also small 
[O] . Typical thermal momenta of single nucleons are pt ~ a/3 m^ Tin . This 



results in corrections of the order [pt/ PiY ~ Tjn/AfTf, about 5% for frag- 
ments like carbon. Final state Coulomb interaction will mainly shift the 
energies upwards, and with it raise the temperature Tf slightly. 

The exact value of Tf as a function of Tin depends on Ep, which in turn 
depends on the freeze-out density as Ep{p) = Ep{p = pQ) (p/po)^''^- It is, per- 
haps, instructive to insert some typical numbers: If we assume a freeze-out 
density of 0.5 po, and a temperature of Tjn = 6 MeV, then we obtain a Fermi 
energy of 24 MeV, and a slope constant (ignoring the recoil correction) of 
Tf = 12 MeV. This number is of the order obtained by projectile fragmen- 
tation reactions, and it therefore is possible that this kind of fragmentation 
reactions are not really a cold breakup of nuclei at nuclear matter density (as 
is conventionally assumed), but the fragmentation of a more dilute system 
at higher internal temperature. 

In summary, by application of the single particle model it was shown 
that it is possible to understand the high values of apparent fragment kinetic 
energy spectra temperatures as compared to temperatures extracted from 
isotope ratios or level population ratios in proton-induced or reverse kine- 
matics fragmentation reactions. One may speculate that this basic mecha- 
nism can also qualitatively explain similar observations for heavy ion induced 
fragmentation events with sizeable compression (by inclusion of radial fiow 
momentum components), but the model employed here might be too simple 
for this case. 
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Figure Captions 

Fig. 1 Fragment kinetic energy spectra for different masses (histograms) 
from solving Eq. ^ witli a value for the variance of 0.014 GeV^. The 
sohd hne is the Gauss limit, Eq. ^. 

Fig. 2 Apparent temperature of fragment kinetic energy spectra (in units 
of the Fermi energy) as a function of the temperature, Tin of the Fermi 
gas. Solid line: numerical solution of Eq. || inserted into Eq. |^. Dashed 
line: analytic approximation, Eq. [TT|. 



